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Abstract 

A fully coupled transient heat and moisture transport in a masonry struc- 
ture is examined in this paper. Supported by several successful applications 
in civil engineering the nonlinear diffusion model proposed by Kiinzel [H is 
adopted in the present study. A strong material heterogeneity together with 
O . a significant dependence of the model parameters on initial conditions as well 

as the gradients of heat and moisture fields vindicates the use of a hierarchi- 
^ ^ cal modeling strategy to solve the problem of this kind. Attention is limited 

1/-^ ■ to the classical first order homogenization in a spatial domain developed 

^ ■ here in the framework of a two step (meso-macro) multi-scale computational 

. scheme (FE^ problem). Several illustrative examples are presented to inves- 

tigate the influence of transient flow at the level of constituents (meso-scale) 
on the macroscopic response including the effect of macro-scale boundary 
conditions. A two-dimensional section of Charles Bridge subjected to actual 
climatic conditions is analyzed next to confirm the suitability of algorithmic 
■ '-3 , format of FE^ scheme for the parallel computing. 
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1. Introduction 

Historical masonry structures all over the world enjoy constant atten- 
tion by many entities including technical audience and public authorities. 
When subject to restoration these structures naturally invite complex anal- 
yses combining both experimental work ^ and numerical simulations jsj. 
Even after closing the scheduled reconstruction steps a continuation of in- 
situ monitoring is now becoming almost standard allowing the engineers not 
only to evaluate the current state of the structure but also to improve the 
predictive capability of theoretical models being supported by up to date 
material data and instantaneous measurements of the state of stress and de- 
formation. A particularly vivid example of this line of inquiry is the Charles 
Bridge information system [^] integrating detailed geometrical description 
of the bridge including changes resulting from previous reconstructions, his- 
torical and contemporary materials forming the bridge as well as novel ma- 
terials and technologies used to improve its durability and serviceable hfe. 
The available results of long-term measurements of temperature and mois- 
ture fields, time-dependent displacement data, complemented by advanced 
methods of computational analysis then open the way to incorporate multi- 
physics, multi-scale, time-dependent and three-dimensional aspects of the 
problem to create a realistic computational model of such a complex struc- 
ture. 

These issues have been supported by our recent study [5i] devoted to the 
homogenization of masonry walls with emphases on the effect of imperfect 
hydraulic contact. On the one hand, it has been demonstrated that the intro- 
duction of interface transition zone at the mesostructural level, considerably 
complicating the computational matter, is essentially negligible for the pre- 
diction of effective properties. On the other hand, an accompanied paramet- 
ric study has shown a substantial dependence of the homogenized properties 
on both the initial and loading conditions. Keeping the format of a simpli- 
fied uncoupled multi-scale scheme with several successful applications partic- 
ularly in civil engineering [sl, [g], 0, to cite a few] would require performing the 
homogenization analysis for various values of water content and certain ref- 
erenced initial values of temperature and relative humidity to construct the 
homogenized macro-scale retention curves. Utilizing these curves in an inde- 
pendent macroscopic study would considerably increase the computational 
efficiency since avoiding a time consuming down scaling for the parameters 
update at every macroscopic time step. Unfortunately, an observed strong 
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dependence of the effective properties on the apphed macroscopic gradients 
may lead to an enormous database of these response functions essentially 
loosing the advantage over a full-fledged coupled multi-scale framework. This 
issue together with the possibility of including the mesostructural morphol- 
ogy and mesostructural material behavior in the macro- level, where typical 
structures are analyzed, without the need for assigning the fine-scale details 
to the entire structure thus motivated the developments in this paper towards 
an iterative FE^ algorithm much similar to that presented in [8| . 

Owing to the finite size of the representative volume element on meso- 
scale we adopt a variationally consistent homogenization developed in jsl to 
reflect a certain size dependency of distributions of macroscopic fields due to 
higher order terms appearing on the left hand side of macroscopic balance 
equations when a non-linear transient flow is assumed on both the macro 
and meso-scale. A short numerical study of this topic is presented in Sec- 
tion 14.11 preceded by the theoretical formulation in Section [31 Section 14.21 
then investigates the influence of macroscopic finite element mesh and appli- 
cation of the associated loading conditions in FE^ scheme on the accuracy 
of evolution of macroscopic moisture and temperature fields. The obtained 
results are then utilized in Section 14.31 devoted to the parallel format of the 
underlying multi-scale analysis performed for a two-dimensional section of 
Charles Bridge. Summary of the essential findings is available in Section |5l 
To keep the paper self-contained a short overview of the adopted constitutive 
model is provided in Section |2J 

In the following text, a and A denote a vector and a symmetric second- 
order tensor, respectively. The symbol V = {d/dx,d/dy,d/dz}'^ stands for 
the gradient representation. All materials are assumed locally isotropic. 



2. Material model 

The literature offers a manifold of material models that allow for the de- 
scription of coupled heat and moisture transport. An extensive overview 
of transport models is presented in the monograph by Cerny and Rov- 
namkova Among others the non-linear diffusion model proposed by 
Kiinzel holds a great potential for an accurate description of transport pro- 
cesses in building engineering including masonry structures and will be adopted 
here to follow up our previous works in this area 10. [sf. 



Kiinzel derived the coupled system of energy and mass balance equations 
based on concepts put forward by Krischer and Kiessl, see e.g. 
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Krischer identified two transport mechanisms for material moisture, one be- 
ing tfie vapor diffusion and the other being described as a capillary water 
movement. In other words, Krischer introduced the gradient of partial pres- 
sure in the air as the driving force for the water vapor transport and the 
gradient of liquid moisture content as the driving force for the water trans- 
port. Kiessl further extended the diffusion model of Krischer and developed 
in 13| his own original version. The unification for the description of mois- 
ture transport in the hygroscopic ip < 0.9 and over hygroscopic ip > 0.9 range 
{(fi is the relative humidity) was achieved with the help of moisture potential, 
which brought several advantages particularly a very simple expression for 
the moisture transport across interfaces. On the other hand, the definition of 
moisture potential in the overhygroscopic range was too artificial, and Kiessl 
introduced it without any theoretical background, see j^. 

For the description of simultaneous water and water vapor transport 
Kiinzel neglected the liquid water and water vapor convection driven by 
gravity and total pressure as well as enthalpy changes due to liquid flow and 
choose the relative humidity as the only moisture potential for both hy- 
groscopic and overhygroscopic range. He also divided overhygroscopic region 
into two sub-ranges - capillary water region and supersaturated region, where 
different conditions for water and water vapor transport are considered. In 
comparison with Kiessl's or Krischer's model Kiinzel's model introduces sev- 
eral simplifications. Nevertheless, the proposed model describes all substan- 
tial phenomena of the heat and moisture transport in building materials and 
the predicted results comply well with the experimentally obtained data [sj. 

Employing the classical Pick's law for the description of water vapor dif- 
fusion, Kelvin's law to simulate the transport of liquid water and Fourier's 
law to account for the flow of heat energy the Kiinzel model integrates into 
the energy balance equation 

"^^"^^ V^[XVe] + KV^[6,V{vpU9)}], (1) 



d^ dt 

and the conservation of mass equation 

= ^"[^-^^] + ^^l^p^iVP^^m] ' (2) 

where H is the enthalpy of the moist building material, w is the water con- 
tent of the building material, A is the thermal conductivity, D^, is the liquid 
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conduction coefficient, Sp is the water vapor permeability, hy is the evapora- 
tion enthalpy of the water, psat is the water vapor saturation pressure, 9 is 
the temperature and ip is the relative humidity. The material parameters, 
both measured and functionally derived, that enter the above equations are 



summarized for the sake of completeness in [Appendix A[ see also ll|, [14 



for more detailed discussion on this subject. Note that the second term on 
the right hand side of Eq. ([T]) represents the change of enthalpy due to phase 
transition being considered the only heat source or sink. 



^i^j 5i> 




nuTfl, r^' ur. 



Figure 1: A two dimensional region 

Application of principal of virtual work then yields the weak form of these 
two balance equations 



d9 dt 



{V69} 



T 



dp., 



sat 



d^ 



V9 



dQ + 



+ [ {V6ip}'^ [K{5pPs.t'Vip}]dn- [ 59q,dT = 0, 
Jn Jrl 



(3) 



^ dw dip 
n dt 



dn + + [ { V(5<^}T [{D^Vif} + {<5pPsatV(^}] dn + 

Jn 

5ipg,dT = 0, (4) 



+ I {V69}'^ i^6py,^V9Un 



p9 

i 1^ 



to be solved numerically for the prescribed boundary and initial conditions, 
see Fig. [T]for the definition of prescribed boundary terms and domain repre- 
sentation. Owing to a strong nonlinear dependence of material parameters 



on both temperature and moisture fields, recall [Appendix A[ the Newton 



Raphson method is generally needed to solve the resulting discretized system 
of equations. 
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3. First order homogenization of non-stationary coupled heat and 
moisture transport 



The present section derives the governing equations of the coupled heat 
and moisture transport in the framework of coupled two-scale analysis of 
FE^ type. In this regard it is presumed that the homogenized macro-scale 
fields are found from the solution of a certain sub-scale (meso-scale) problem 
performed on a representative volume element (RVE) being identical, at least 
in a statistical sense, to a real meso-structure from both the geometrical and 
material composition point of view. Such an RVE is then usually termed 
the statistically equivalent periodic unit cell (SEPUC) [Hj. Examples of 
such SEPUCs for both regular as well as irregular masonry walls adopted 
herein are plotted in Fig. |6]^d)(e). It has been advocated in [8] that for 
a finite size RVE the assumption of transient fiow on both the macro and 
meso-scale introduces certain non-local, size dependent, terms in equations 
governing the macroscopic response. Some numerical simulations addressing 
this issue are presented in Section 14. whilst the theoretical grounds are 
provided next following closely, although in more abbreviated format, the 
variationally consistent homogenization outlined in detail in jsj. 

To introduce this subject suppose that a local field a can be replaced by 
a spatially homogenized one (a) such that 

/ adn ^ [ {a)^dn = [ (-^ [ adnn]dn, (5) 
Jn Jn Jn V Jn^ J 

j^adT ^ ^(a)^dr = ^(^^^ adFn^dr, (6) 

where Qn and Fn represent the internal and boundary parts of the SEPUC. 
In what follows, owing to the space limitation, we shall treat only the energy 
balance equation (|3]) which upon employing Eqs. ([5]) and ([6]) becomes 



dn 

n 



+ [ ({vM"^[/^45pPsatVv9}]\ d^]- /" {6eq,)^dT = o. 

Jn^ Jp? 

In the spirit of the first order homogenization it is assumed that the 
macroscopic temperature and relative humidity vary only linearly over the 
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SEPUC. This can be achieved by loading its boundary by the prescribed tem- 
perature 0^°™ and relative humidity ^^'^"^ derived from the uniform macro- 
scopic temperature VB and relative humidity V*!* gradients. In such a case 
the local temperature and relative humidity inside the SEPUC admit the 
following decomposition 



e{x) = Q{x^) + {VQ}^{x-x^} + e*{x) 

^{x) = $(X°) + {V<l>}"^{a;-X°} + (^*(a;) 



Q^°'^{x) + e*{x), 

-- <^^°'^{x) + ^*{x) 



(8) 
(9) 



where 9*{x) and f*{x) are the fluctuations of local fields superimposed onto 
linearly varying quantities 6^°™ (a;) and $^°'^(a;) . The temperature e(X°) 
and the moisture at the reference point are introduced to link the 

local fields to their macroscopic counterparts. For convenience the SEPUC is 
typically centered at X^. Henceforth, the local fluctuations will be demanded 
to be periodic, i.e. the same values are enforced on the opposite sides of a 
rectangular SEPUC. 

Next, substituting Eqs. ([8]) and ([9]) into Eq. ([7]) and collecting the terms 
corresponding to 50*^°™, 5$'^°™ and 66*, 6(p* splits the original problem ([3]) 
into the homogenized (macro-scale) problem 



dH d9 



n \ 



{W9} + K 



+ 



{V5($^°-)}"^[/^.{^pPsatV^}]\ dn- f ({5e'^n^.>ndr = o, 



and the local sub-scale (meso-scale) problem 



de 



dQ + 



(10) 



+ 



{{69*}q^)^dT 



dQ + 



=0 duo to poriodicity 



{69*}q^d{dna)dn = 0, 



which is satisfied identically owing to the assumed periodic boundary condi- 
tions. Solving Eq. ( ITTj) for the prescribed increments of V0 and V$ provides 
the instantaneous effective properties and storage terms that appear in the 



7 



macro-scale equation ( ITU]) . Because of a strong non-linearity the two equa- 
tions must be solved iteratively in a certain nested loop, see 16, sj for further 
reference. 

Since details on the solution of Eq. f|TT]) are available in our preceding 
paper jsf we limit our attention to the macro-scale problem and write the 
first term of Eq. f lTU]) with the help of Eqs. (|H]) and as 



{50hom| 



dH d9 



{S{e + {VQ}^{x-x'})} 



dH d9 

dedt 



dQ 



{6VQ}'^ {x - X"} 



dH d9 
d^dt 



dfi. 



thus clearly identifying the solution dependence on the actual size of the 
SEPUC through the second term in the integral (fT2|) . We may now substitute 
from Eq. into Eq. to get 



dn- 



c 

{SVQ} 







— s/ 



dt 

T 



{\V9} + K <! ^PV^'^O 



dn 



^ ({5V$}T [K{6pP,,,V^}]^^ dQ + j^^ ({5e}T q}j^ dV = 0. 



(13) 



9cxt 



An analogous approach can be applied also to the moisture transport equa- 
tion (jl]) to arrive, after classical finite element discretization, into a dis- 
cretized system of coupled macroscopic heat and moisture equations 



drfl 



dt 



9e 



(14) 
(15) 



which have to be properly integrated in the time domain adopting for example 
the Crank-Nicolson integration scheme. Details on the numerical implemen- 



tation are available in 14 



(12) 
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4. Examples 



Several illustrative example problems were analyzed to address the non- 
linear transient coupled heat and moisture transport assumed on both scales, 
the influence of the way of prescribing the macroscopic loading conditions 
closely related to the macro-scale flnite element mesh and flnally the solution 
strategy exploiting the parallel computation. The same material data were 
adopted in all analyses. These were obtained from a set of experimental 
measurements providing the hygric and thermal properties of mortars and 
bricks/stones, which have been used in the reconstructions works of historical 



buildings in the Czech Republic including Charles Bridge, see [17|. The 
measured material parameters of individual masonry phases listed in Table. [T] 
then served to derive the non-measurable transport coefficients presented in 
Eqs. (II]) and ([2]), see also Appendix A 



parameter 






brick 


mortar 


Wf [kgm" 




free water saturation 


229.30 


160.00 


wso [kgm- 




water content at ip = 0.8 [-] 


141.68 


22.72 


Ao [Wm 




thermal conductivity 


0.25 


0.45 


fetes [-] 




thermal conductivity supplement 


10 


9 


ps [kgm" 




bulk density 


1690 


1670 


[-] 




water vapor diffusion resistance 


16.80 


9.63 


A [kgm" 


-2g-0.5] 


water absorption coefficient 


0.51 


0.82 


Cs [3kg- 




specific heat capacity 


840 


1000 



Table 1: Material parameters of individual phases. 



4.1. Influence of transient flow at meso-level 

This section supports through numerical simulations the theoretically pre- 
dicted size dependence of the homogenized response first suggested for the 
case of a non-linear single variable diffusion problem in [sj and also estab- 
lished here in Section [3] for the special case of the coupled heat and moisture 
problem in the framework of Kiinzel's constitutive model. 

In doing so we considered three particular units cells in Fig. [2](a) varying 
in size from millimeters to decimeters. Each cell was loaded by the same 
constant gradients of temperature and moisture along the x-direction. The 
resulting evolutions of the fluctuation part of the local temperature at the 
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center of individual cells appear in Fig. I21^b) clearly manifesting the influ- 
ence of the size of the cell which necessary projects into the prediction of 
the homogenized properties and thus evolution of the predicted macroscopic 
response. For the largest cell the steady state was reached in about 80 [h]. 
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Figure 2: (a) Investigated periodic unit cells, (b) resulting evolution of fluctuation tem- 
perature 



4-2. Influence of macro structural finite element mesh 

One of the concerns of the implementation of FE^ scheme is the appli- 
cation of boundary conditions on the macro-scale keeping in mind the scale 
transition requirement and periodic boundary conditions imposed on the 
meso-scale. This may become important particularly with a relatively large 
periodic unit cells which may even exceed the size of macroscopic elements 
in the vicinity of outer boundary where the macro-elements should be fine 
enough to ensure a smooth and accurate evolution of driving variables from 
the outside into the inner parts of a structure. 

To address this issue we studied two types of macro-scale discretizations 
adopted in the two-scale (meso-macro) analysis. The corresponding macro- 
scale finite element meshes appear in Figs. [3]^a),(b). The case when the fine- 
scale details are assigned to the entire structure is plotted in Fig. El^c). This 
mesh, consisting of 7050 triangular finite elements, served to evaluate the ac- 
curacy of the two former discretizations. Fig. |3]^a) shows 108 macro-elements 
each representing a single meso-problem with assigned periodic boundary 
conditions, whereas the case in Fig. [3t^b) assumes the outer boundary being 
fully discretized (note that only one-directional flow is considered). There, 
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only the inner part consisting of 72 elements is subject to multi-scale analysis 
whilst the outer part is modeled as a structure with a real masonry bonding 
consisting of 3888 finite elements. A multi-point constrains were introduced 
to account for an incompatible discretization along the common interface. 
This latter case is, therefore, expected to heal inaccuracies in the estimation 
of temperature and moisture fields in the region close to the surface layer. 























^ -H ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 






















^ —H ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 

1 ' 1 !■ I 1 1 1 1 1 1 1 ' 1 






















^ 1 ' '■■[ ■'■[ 'l 'l '■ 
























1 ' 1 1 1 1 1 1 1 1 1 
























^ -T-' ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 






















^ —H ' I ' I ' 1 ' 1 ' 1 ' I 












































! ' ! ' ; ' ; ' ; ' 1 ' 



(a) (b) (c) 



Figure 3: Different finite element representations of masonry wall {Ix — 1.92 [m], ly = 
1.80 [m]) - (a) full multi-scale scheme, (b) semi multi-scale scheme, (c) full fine-scale dis- 
cretization 

The following boundary conditions were imposed: on the right-hand side 
(interior) a constant temperature of 24 [°C] and a constant relative humid- 
ity 0.5 [— ] were maintained, while on the left-hand side (exterior) the real 
climatic data collected over the entire year were prescribed, see Fig. |H 




Figure 4: Annual loading conditions - (a) temperature, (b) moisture 
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The results appear in Fig. [S] showing variation of the temperature and 
moisture along the mid section of the wall after the duration of load of 10 [h] 
and 100 [h], respectively, derived for the macroscopic time step equal to 1 [h]. 
Clearly, the notable difference between the exact (full fine-scale discretiza- 
tion) and full multi-scale scheme can be observed in surface layers only and 
this difference almost disappears with a sufficiently long duration of time. It 
thus appears that the refined representation of the surface layer through the 
semi multi-scale scheme, although more accurate compare to full multi-scale 
scheme, does not bring any particular advantage. This is supported by the 
calculated average and absolute errors stored in Table |2] taking into account 
all nodal macroscopic temperatures and moistures in the domain over all 
time integration steps. Note that for the sake of comparison the fine-scale 
variables (solution employing the mesh in Fig. EJj^c)) were averaged over the 
cell basically covered by two macro-elements in Fig. |3]^a). 




(d) 



Figure 5: Comparison of different niacrostructural computations - (a) temperature profile 
after t = 10 [h], (b) moisture profile after t = 10 [h] (c) temperature profile after t — 100 [h], 
(d) moisture profile after t — 100 [h] 

The presented results thus promote the more accurate semi multi-scale 
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scheme only for calculations demanding higher accuracy of local results es- 
pecially in initial stages of computation and/or examples with fast changing 
boundary conditions. 



type of comparison 




calculation of temperature 
- fine-scale vs. multi-scale 



2.62 
0.71 



0.11 
0.03 



- fine-scale vs. semi multi-scale 
calculation of relative humidity 

- fine-scale vs. multi-scale 



0.18 
0.05 



0.001 
0.001 



- fine- scale vs. semi mult i- scale 



Table 2: Averaged relative and absolute errors 



4-3. Parallel computation 

The essential request by the contractor when studying the mechanical re- 
sponse of Charles Bridge to provide the basis for reconstruction works was a 
full scale three-dimensional analysis of the bridge. Performing such an anal- 
ysis in a fully coupled format on a single computer would be computationally 
unfeasible thus creating the need for a parallel computing. Concentrating on 
the implementation part of the parallel version of FE^ scheme we limit our 
attention to a two-dimensional section of Charles Bridge subjected, however, 
to real climatic data displayed already in Fig. m Extension to a fully three- 
dimensional problem is under current investigation and will be presented 
elsewhere. 

As already discussed in the previous section, the present FE^ based multi- 
scale analysis assumes each macroscopic integration point be connected with 
a certain mesoscopic problem represented by an appropriate periodic unit 
cell. The solution of a meso-scale problem then provides instantaneous effec- 
tive data needed on the macro-scale. Such an analysis is particularly suitable 
for a parallel computing because the amount of transferred data is small. In 
this regard, the master-slave strategy can be efficiently exploited. To that 
end, the macro-problem is assigned to the master processor while the solu- 
tion at the meso-level is carried out on slave processors. At each time step 
the current temperature and moisture together with the increments of their 
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gradients at a given macroscopic integration point are passed to the slave pro- 
cessor (imposed onto the associated periodic cell), which, upon completing 
the small scale analysis, sends the homogenized data (effective conductivities, 
averaged storage terms and ffuxes) back to the master processor. 

If the meso-scale problems are large enough, the ideal solution is to as- 
sign one meso-problem to one slave processor. Clearly, even for very small 
macro-problems with a few thousands of finite elements, the hardware re- 
quirements would be in such a case excessive. On the other hand, if the 
meso-problems are relatively small, i.e. they contain small number of finite 
elements, the corresponding analysis might be even shorter than the data 
transfer between the processors. Then, the computational time associated 
with the data transfer between the master processor and many slave proces- 
sors may grow excessively. It is worth mentioning that this time consists of 
two contributions. The first one represents the latency time (the processors 
make connection) which is independent of the amount of transferred data 
whilst the second contribution clearly depends on the amount of data being 
transferred. For small meso-problems it is therefore reasonable to assign sev- 
eral of them to a single slave processor. The master processor then sends a 
larger package of data from many macroscopic integration points at the same 
time to each slave processor so that the latency time does not play a crucial 
role. This approach was adopted hereinafter. 

Fig. El^a) displays a three-dimensional segment of Charles Bridge exam- 
ined in the original three-dimensional static calculation A two-dimensional 
cut through the mid part of a Charles Bridge arch examined for the paral- 
lel computation appears in Fig. Mjo) together with four material regions. 
These are associated with two heterogeneity systems of the bridge, one rep- 
resenting a regular sand stone masonry of side walls, fence and arches and the 
other corresponding to an irregular quarry masonry made of arenaceous marl 
blocks and sand and black lime mortar filling the inner part of the bridge. 
For simplicity, the bridge deck was assigned the regular pattern. The corre- 
sponding periodic unit cells employed for the meso-scale analysis are plotted 
in Figs. El^d) andEl^e), respectively. 

The finite element mesh used at the macro-level is evident from Fig. M^c) 
featuring 7, 081 nodes and 13, 794 triangular elements with a single integra- 
tion point thus amounting to the solution of 13, 794 meso-problems at each 
macroscopic time step. This figure also shows decomposition of the macro- 
problem into 12 slave processors. The numbers of elements in individual 
sub-domains being equal to the number of meso-problems handled by the 
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Figure 6: (a) Three-dimensional view of Charles Bridge with a two-dimensional A-A 
section analyzed, (b) Analyzed section showing material regions with assigned meso-scale 
unit cells, (c) Macrostructural mesh with identified loading conditions and decomposition 
into sub-domains representing individual slave processors (Lj: — 10.40 [m], ly = 6.82 [m]), 
(d) Mesostructural mesh of regular bonding of masonry (SEPUC 1: Ix — 0.45 [m], ly = 
0.44 [m]), (e) Mesostructural mesh of irregular quarry masonry (SEPUC 2: = 0.45 [m], 
Z;, = 0.35 [m]). 



assigned slave processor are listed in Table [31 It should be noted that the as- 
sumed decomposition of the macro-problem is not ideal. In comparison with 
domain decomposition methods, the macro-problem has to be split with re- 
spect to the heterogeneity of the material resulting in the variation of number 
of elements in individual sub-domains between 1046 and 1748. 

The number of elements in the two meso-problems amounts to 265 (160 
nodes) for SEPUC 1 and to 414 (239 nodes) for SEPUC 2, respectively. 
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Processor No. 


1 


2 


3 


4 


5 


6 


7 8 




9 


10 


11 


12 








No. of 


1218 


1748 


1046 


1052 


1214 


1210 


1052 1054 


meso-problems 


1046 


1052 


1054 


1048 









Table 3: Decomposition of the macro-problem into sub-domains. 



Similarly to the macro-problem, the meso-problems have to account for the 
material heterogeneity. Clearly, the ideal speedup and load balancing are 
obtained when the decomposition of the macro-problem reflects the meso- 
problem meshes. However, this is considerably more difficult when compared 
to the classical mesh decomposition. 

The actual analysis was performed on a cluster built at our department. 
Each node of the cluster is a single processor personal computer Dell Optiplex 
GX620 equipped with 3.54 GB of RAM. The processors are Intel Pentium 
with the frequency 3.4 GHz. The cluster is based on Debian linux 5.0 and 
32-bit architecture. 




Figure 7: Evolution of macroscopic (a) temperature, (b) moisture 

Although various mesoscopic heterogeneity patterns were properly ac- 
counted for, the material data of individual constituents (bricks or stones 
and mortar) were taken the same for both the outer and inner part of the 
bridge, recall Table [T]for specific values. The initial conditions on the macro- 
scale were set equal to (p = 0.5 [-] and 6^ = 20 [°C] and the year round 
variation of moisture and temperature in Fig. H] was imposed onto all outer 
surfaces of the bridge, see Fig. MyC). Owing to the computational demands 
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the macroscopic time increment was set equal to 10 [h], which agreed with 
the real computational time equal to 2.08 minute for each time step. In 
view of the results presented in Section 14. 2[ recall Fig. [5l this justifies, al- 
though at the loss of accuracy at the initial stage of computation, the use of 
full multi-scale scheme. One particular example of the resulting evolutions 
of macroscopic temperature and moisture at the selected nodes labeled in 
Fig. ^c) is seen in Fig. [71 

5. Conclusions 

A fully coupled multi-scale analysis of simultaneous heat and moisture 
transport in masonry structures was implemented in the framework of FE^ 
computational strategy. Two particular issues were addressed: the influence 
of the finite size of SEPUC when running the transient analysis on both scales 
and the way of introducing loading on the macro-scale. While the former one 
plays a significant role in the estimates of macroscopic response, the latter 
one proves important only in the initial stages of loading. 

Special attention was accorded to the implementation of FE^ concept 
in the parallel format employing the master-slave approach. Although not 
qualitatively fully acceptable, the two-dimensional example of Charles Bridge 
raised a number of questions to the solution efficiency particularly with ref- 
erence to a proper subdivision of the analyzed macro-domain and local finite 
element mesh of individual meso-scale SEPUCs. The present findings sum- 
marized in Section 14.31 will be utilized in a fully three-dimensional analysis 
being the subject of our current research effort. 
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Appendix A. 

The list of material parameters to be obtained experimentally are stored 
in Tabled] The transport coefficients that enter Eqs. ([T]) and ([2]) are provided 
by 
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w - water content [kgm 

"-Z^' (A.1) 

where wj is the free water saturation and h is the approximation fac- 
tor, which must always be greater than one. It can be determined 
from the equihbrium water content (wgo) at 0.8 [-] relative humidity by 
substituting the corresponding numerical values in equation flA.ip . 

5p - water vapor permeability [kgm~^s^^Pa^^], 

5, = (A.2) 

where n is the water vapor diffusion resistance factor and 6 is the vapor 
diffusion coefficient in air [kgm~^s~^Pa~^] given by 

2.306 -10-^ p. / g + 273.15 y-'^ 
{6 + 273.15) p V 273.15 J ' ^ ' 

with p set equal to atmospheric pressure pa = 101325 [Pa] and = 
R/Mu, = 461.5 [Jkg-^K-^]; i?isthe gas constant (8314.41 [Jmor^K-^]) 
and My, is the molar mass of water (18.01528 [kgmol^^]). 



= D^ — , (A.4) 



- liquid conduction coefficient [kgm ^s ^ 

dw 

where Dy^ is the capillary transport coefficient given by 

D^ = 3,s(—\ .io3WK-i)^ (A.5) 

where the derivative of the moisture storage function is obtained 

dip 

by differentiating Eq. (lA.ll) . 

A - thermal conductivity [Wm~^K~^], 

A = Ao (l + ^) , (A.6) 

where Aq is the thermal conductivity of dry building material, pg is the 
bulk density and 6tcs is the thermal conductivity supplement. 
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Psat - water vapor saturation pressure [Pa], 



Psat = 611 exp ( I ' ^^■^'> 



where 



(A.8) 



a = 22.44 00 = 272.44 [°C] < [°C] 

a = 17.08 00 = 234.18 [°C] e>0 [°C] 
• hy - evaporation enthalpy of water [Jkg~^] 

/273 1 ^\ (0-167+3. 6710-*6») 

/i, = 2.5008 ■ lOM —y- j . (A.9) 
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